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ABSTRACT 


It  has  become  of  recent  interest  to  complement  electromagnetic  responses  in  intensity  by 
the  additional  consideration  of  their  polarizations.  In  order  to  affect  the  latter,  media  have  been 
considered  that  contain  inclusions  of  randomly  distributed  and  oriented  short,  thin  wires  (needles), 
whose  polarized  scattering  amplitudes  can  modify  the  well-known  polarization  of  surface 
reflections  that  is  governed  simply  by  the  so-called  Fresnel  coefficients.  In  the  present  study,  we 
analytically  obtain,  for  the  case  of  a  multiplicity  of  perfectly  conducting  wires  imbedded  in  a 
dielectric,  the  corresponding  wave  scattering  amplitudes  and  their  superposition  with  the 
specularly  surface-reflected  wave  amplitude  when  both  are  being  generated  by  the  same  incident 
wave.  This  allows  a  calculation  of  the  total  retumed-wave  polarizations  which  we  express  by  the 
real  Stokes  parameters. 


ADMINISTRATIVE  INFORMATION 

Professor  Oberall  conducted  this  investigation  under  the  auspices  of  the  American  Society 
of  Engineering  Educators  (ASEE).  The  ASEE  program  at  the  Carderock  Division,  Naval  Surface 
Warfare  Center  is  administered  by  F.  Halsall,  Code  01 12.  The  purpose  of  the  continuing  effort  is 
to  study  the  fundamental  interactions  of  materials  and  electromagnetic  energy. 
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INTRODUCTION 


The  reflected  return  of  electromagnetic  waves  fi’om  dielectric^’^  (including  lossy  ones)  or 
metallic  surfaces^  is  well-known,  and  can  be  analytically  expressed  by  the  Fresnel  coefficients*’'* 
which  also  describe  the  refi-acted  wave  that  penetrates  into  the  medium.  These  expressions  are 
valid  for  infinite  homogeneous,  isotropic  media  with  plane  surfaces  that  are  smooth  compared  to 
the  wavelength.  These  expressions  also  describe  the  polarization  of  the  reflected  and  refracted 
waves.  If  the  polarization  (electric)  vector  of  the  incident  wave  is  perpendicular  to  the  plane  of 
incidence  (this  case  termed  “r”),  the  same  will  hold  for  the  reflected  and  refracted  waves.  For  the 
incident  polarization  vector  that  is  parallel  to  the  plane  of  incidence  (case  termed  ‘7’)»  reflected 
and  refracted  wave  polarization  vector  will  be  parallel  to  the  plane  of  incidence.  Thus  no  coupling 
of  linear  polarizations  occurs.  For  the  ^case,  the  reflected  wave  vanishes  at  the  Brewster  angle  of 
incidence,  so  that  for  a  plane  incident  wave  of  mixed  polarization,  the  reflected  wave  is  of  purely 
r-type  at  that  angle. 

One  may  attempt  to  modify  the  polarization  of  waves  reflected  at  the  surface  of  a 
dielectric  (which  may  be  lossy),  by  imbedding  short,  thin  wires  (needles)  in  it  that  are  randomly 
distributed  and  oriented.  Our  previous  study*  has  considered  such  media  for  the  case  of 
inclusions  that  are  small  compared  to  the  wavelength;  for  that  case,  one  was  able  to  use  “effective 
medium  theory”®  which  represents  the  material  as  one  homogeneous  medium  with  effective 
properties.  In  the  present  case,  we  chose  to  consider  wires  long  compared  to  the  wavelength  (and 
also  long  compared  to  their  diameter),  so  that  effective  medium  theory  is  no  longer  applicable. 

The  corresponding  scattering  cross  section  of  one  wire  can  be  obtained  fi’om  Reference  4,  and  the 
scattered  wave  combines  coherently  with  the  surface-reflected  wave  (given  by  the  corresponding 
Fresnel  coefBcient),  both  being  generated  by  the  same  incident  wave.  The  effects  of  randomly 


2 


distributed  multiple  wires  can,  however,  be  added  together  incoherently  as  shown  by  van  de 
Hulst’.  Polarization  of  the  scattered  wave  is  purely  of  ^type  and  thus  adds  to  the  ^type  Fresnel 
reflection  while  the  r-type  Fresnel  reflection  is  not  affected  by  the  wire  scattering.  The  total  ^-type 
response  (reflected  and  scattered)  is  found  to  be  superposed  with  different  phases,  rendering  its 
amplitude  complex  so  that  the  return  polarization  cannot  be  advantageously  described  by  the 
customary  simple  polarization  formula  (difference  of  r  and  1  amplitudes  divided  by  their  sum), 
which  would  now  become  complex.  Instead,  we  characterize  polarizations  by  the  real 
components  of  the  Stokes  4-vector,’  of  which  only  the  first  two  are  non-vanishing  since  only 
linear  polarizations  are  involved  here.  Numerical  calculations  show  the  influence  of  scattering 
fi'om  the  imbedded  wires  on  the  polarization  of  the  total  (reflected  and  scattered)  returns. 

REFLECTED  AND  REFRACTED  WAVES 

We  consider  a  half  space  (which  later  on  will  be  assumed  to  have  a  finite-size  surface  area 
A)  filled  with  a  (possibly  lossy)  dielectric  that  contains  a  total  of  N  randomly  distributed, 
randomly  oriented,  perfectly  conducting  cylindrical  needles  of  length  l=2h  and  diameter  2a,  with 
a«h.  (Later  on,  these  needles  can  be  taken  as  being  of  fiunite  conductivity  and  permittivity,  which 
may  be  the  subject  of  a  subsequent  analysis).  In  practice,  the  half  space  will  have  a  finite  depth, 
and  N  is  the  number  of  needles  contained  in  its  volume  (A  x  depth). 

The  physical  situation  is  as  shown  in  Fig.  1.  The  wave  E  =  Eo  exp  ik-r  is  incident  from 
below  (at  angle  Si)  through  a  medium  with  permittivity  e  and  permeability  p.  (we  shall  always 
take  the  case  p  =  1),  considered  to  be  air.  The  dielectric  half  space  (top)  has  permittivity  e'  and 
permeability  p'  (again  p'=l)  and  contains  an  upgoing  refracted  wave  E'  (at  angle  Sr)  while  also 
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generating  a  downgoing  reflected  wave  E"  (at  angle  dj).  The  various  waves,  electric  vectors  E 
and  propagation  vectors  k  are  indicated  in  the  figure. 

The  needle  scattering  of  the  refi'acted  wave  E',  which  is  also  indicated  in  Fig.  1,  will  be 
discussed  in  the  following  section.  At  this  time,  we  shall  complete  the  discussion  of  the  refi’acted 
and  reflected  waves  shown  in  Fig.  1 .  The  polarization  of  the  incident  wave  is  assumed  to  be 
purely  linear,  and  we  shall  distinguish  the  two  cases  Er  and  E,as  indicated  in  Fig.  1 .  These 
amplitudes  are  designated  as  Eo"  or  Eo^  while  the  phase  factor  of  the  incident  wave  is  exp  (ik  r). 
As  mentioned  before,  an  incident  polarization  E,  generates  Ev  and  E",  only,  and  an  incident  E^ 
generates  E',  and  E"<  only,  see  Ref  1 .  For  the  refi’acted  amplitudes  Eo'  and  reflected  amplitudes 
Eo",  one  has^  as  expressed  in  terms  of  the  Fresnel  coefficients  F„fr  and  Fka,  respectively: 


*  refr  ^0 

(la) 

F" 

•■^0  ^  rcfl  ^0 

(lb) 

The  indices  of  refraction  are  n  =  and  n ' 

=  (p'e')*'^,  and  the  propagation  vectors  have  the 

magnitudes 

1  k  1  =  1  k"  1  =  k  =  (©/c)  (pe)'^  =  n©/c 

(2a) 

1  k'  1  =  k'  =  (©/c)  (p'e')^^  =  n'©/c 

(2b) 

where  ©  is  the  angular  frequency  and  c  the  light  velocity  m  vacuum.  The  angles  of  incidence  di 
and  of  refraction  S,  are  related  by  Snell’s  law, 

sin  Si  /  sin  S,  =  k'/k  =  n'/n.  (3) 

In  terms  of  these  quantities,  the  Fresnel  coefficients  are 
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2ncos&; 


FU= 


F«fl= 


n  cosdit— Vn'"-  n"  sin^dj 

n  cosS-'^  n'^-  n‘  sin^^i 

n  cosdi^-^  Vn'^-n^  sin^di 
M* 


(4a) 


(4b) 


w _ 2nn'  cos  &i 

n'^  cosdjif  nVn'^-n^  sin^Sj 
M* 

^  cosdj-  nV  n'^-n'  sin^S: 

=  - - 

n'^  cosSiH-n-v/n'^-  n'  sin^S: 


(4c) 


(4d) 


For  lossless  media,  all  the  parameters  in  these  expressions  are  real  and  there  is  no  problem 
employing  the  conventional  formula  for  the  polarization  of  the  reflected  wave, 

P„fi  =  (E" Em  Eo" '+  E?0  (5a) 

=  (FU-Flcfi)/(Fkfl  +  FU)  (5b) 

which  in  this  case  is  a  real  quantity.  For  lossy  media,  however,  this  expression  for  Pren  would 
become  complex  and  in  that  case  it  is  advantageous  to  describe  polarizations  by  the  (real)  Stokes 
parameter,  see  the  following  section.  Also,  if  needle  scattering  is  added,  the  amplitudes  in  Eq. 

(5  a)  would  become  complex  even  for  lossless  media,  so  that  here  a  use  of  the  Stokes  parameters 
is  indicated,  as  will  be  shown  in  the  next  section  below. 

SCATTERING  FROM  IMBEDDED  WIRES 
The  scattering  process  of  the  refi’acted  wave  E'  when  incident  on  a  thin  wire  needle 
imbedded  in  the  dielectric  is  indicated  in  Fig.  1 .  The  first  thing  to  be  noted  is  that  a  very  special 
scattering  geometry  is  called  for  in  the  present  situation.  Since  the  effect  of  wire  scattering  on  the 
polarization  of  the  reflected  wave  E"  is  to  be  investigated  here,  the  refi'acted  wave  that  after  the 
scattering  becomes  E'„  and  then  re-emerges  from  the  medium  as  E"  „  must  be  such  that  this 
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emerging  wave  E"c  is  directed  parallel  to  the  reflected  wave  E"  (with  k"  =  k"),  so  that  both  of 
these  (coherently  interfering)  waves  can  be  observed  together,  and  their  combined  polarization 
obtained.  This  forces  us  to  consider  the  symmetric  scattering  geometry  of  Fig.  1  only,  with 
scattered  and  re-emergence  angles  S,  and  Sias  indicated,  and  the  whole  process  taking  place  in 
the  (xz)  plane  of  incidence  exclusively. 

The  geometry  of  scattering  from  a  thin  cylinder  is  defined  in  Reference  4,  and  is  shown 
here  in  Figs.  2(a)  and  (b)  in  a  coordinate  system  (xyz).  Incident  and  scattered  waves  are  assumed 
to  be  directed  at  angles  'Pj  and  from  the  cylinder  normal,  and  the  (linear)  polarizations  of  both 
are  here  kept  more  general  than  those  considered  earlier  in  this  report  (r  or  t),  i.e.,  lying  at  angles 
Yi  (y«)  froin  the  plane  of  incidence  (scattering).  Our  basic  coordinate  system  (xyz),  however,  is  as 
shown  in  Fig.  1,  so  that  if  the  wire  is  directed  at  spherical  angles  (d,  cp)  in  that  system,  a 
coordinate  transformation  between  the  two  systems  is  in  order.  This  is  sketched  in  Fig.  3(a)  for 
the  wave  incident  on  the  wire,  and  in  Fig.  3(b)  for  the  scattered  wave.  The  spherical-triangle 
relationships  then  give  the  equations 

sin  Yi  =  cos  fr  cos  d,  +  sin  S  sin  fr*  cos  (p  (6a) 

sin  Ys  =  cos  fr  cos  frr  -  sin  fr  sin  frr  cos  (p  (6b) 

which  will  be  needed  later  in  this  calculation. 

The  amplitude  of  the  refracted  wave  E'  as  scattered  by  the  wire  needle  imbedded  in  the 

dielectric  (Fig.  1),  may  be  obtained  by  using  the  procedures  outlined  in  the  Radar  Cross  Section 

Handbook.'*  This  reference  only  quotes  the  scattering  cross-section  (squared  amplitude),  but  the 

amplitude  itself  (including  phase)  can  be  found  in  analogy  to  van  de  Hulst’s  derivations^  that  are 

based  on  Huygens’  principle.  One  should  note  that  this  reference  employs  waves 

6 


exp(-ikz+i(ot),  so  that  the  complex  conjugate  of  the  expressions  quoted  therein  must  be  used  for 
our  purposes. 

We  will  need  the  far-field  scattering  amplitude  of  the  wire  which  at  distance  r  behaves  as  a 

0 

spherically  spreading  wave  e  (This  is  obtained  here  schematically  since  the  subsequent 

emergence  of  the  scattered  wave  from  the  medium  will  transform  it  into  a  spherical  wave  f  ^e  ^ 

thereafter.)  This  far-field  amplitude  is  found  by  van  de  Hulst  starting  from  the  near-field 

0 

amplitude  which  spreads  cylindrically  OC  and  transforming  it  into  the  far-field  amplitude 

by  invoking  Huygens’  principle,  which  obtains  the  propagated  field  as  a  superposition  of  spherical 
wavelets  that  emanate  from  each  point  of  the  field  before  its  propagation.  Van  de  Hulst’ s  result 
of  this  approach  is  given  in  Ref  7,  p.  305  as  the  amplitude  u, 

u  =  (2ih/7ir)ET  exp(ik'r)  (7) 

where  E  is  a  strongly-peaked  x*^  sin  x  -  type  function  with  side  lobes,  to  be  defined  later,  and  T  is 
the  near-field  scattering  amplitude  of  the  cylinder  (which  in  this  case  may  be  taken  as  infinitely 
long).  The  choice  of  this  amplitude  is  governed  by  the  relation  between  cylinder  radius  and 
wavelength,  and  we  here  choose  the  amplitude  for  T  given  as  a  closed-formed  expression  by  Ruck 
et  al.  (Ref  4,  p.  21 1),  valid  for  very  thin  perfectly  conducting  cylinders  (k '  a«  1)  in  the  form 
E|f  “  =  E|f(7i/2k'r)^^  [ln(0.8905  k'a)-i7c/2]'*  exp[i(k'r  +  7i/4)]  (8) 

where  the  polarization  is  accounted  for:  E||  are  the  field  components  parallel  to  the  cylinder  axis. 
[Although  Ruck  refers  to  this  expression  as  a  “far  field”,  it  is  here  the  cylindrically  spreading  field 
relatively  close  to  the  cylinder  which  in  van  de  Hulst’ s  approach  is  taken  as  the  “near  field”.]  The 
scattering  amplitude  resulting  from  the  combination  of  Eqs.  (7)  and  (8)  is  not  given  in  Ref  7, 
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where  only  the  bistatic  scattering  cross  section  Ob  is  quoted.  The  latter  is  defined  (by  Ref.  7  and 
the  general  “engineering  literature”)  as 

Ob  =  4n  do/dQ  (9a) 

which  is  written  here  in  terms  of  the  expression  do/dQ  used  in  the  “scientific  literature”: 

do/dQ  =  r^|E«,|7|Ei„cP  (9b) 

The  expression  for  Ob  may  be  calculated  fi'om  Eqs.  (7)  and  (8)  and  compared  with  that  of  Ref  7, 
in  order  to  verify  (up  to  possible  phases)  the  correctness  of  our  calculation.  Inserting  in  Eq.  (9a) 
and  noting  that  (cf  Fig.  2)  projecting  E'  (or  E'“)  on  the  xz  plane  (or  the  scattering  plane) 
introduces  factors  cos  Yi(or  cos  y,)  while  the  further  projection  on  the  z  axis  adds  factors  with  'Fi 
(or  ^,),  the  latter  in  the  manner  E'“  cos  =  E'  cos  Ti  (...)  in  Eq.  (8).  This  leads  to  the  result 
Ob  =  47ch^(cos  Yj/  cos  ^s)^{sin[k'h(sin  'Fi  +sin  %)]/  k'h(sin  'Fj  +sin  Ta)}^ 

•  cos^Yi  cos^Y»/  [ln70.8905  k'a)  +  7t74],  (10) 

in  agreement  with  Ref  4,  p.  304.  Note  that  from  Eq.  (6),  we  have 

sin  Yi  +  sin  'F,  =  2  cos  d  cos  •  (1 1) 

It  should  be  added  that  in  Ref  4,  Eq.  (4. 3-4.4)  our  factor  (cos  Yi  /  cos  'F,)^  appears  upside-down, 
but  we  believe  that  our  version  is  the  correct  one,  as  can  be  backed  up  by  certain  physical 
arguments.  However,  in  our  problem  it  will  be  shown  below  that  Y,=  -T^i,  so  that  the  factor  in 
question  equals  unity  in  either  version. 

The  peaked  angular  function  E  is  now  identified,  using  Eq.  (1 1),  as 
E  =  sin  (2k'h  cos  S  cos  dr)  /  (2k'h  cos  d  cos  dr)  (12) 
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With  a  =  2k'h  cos  dr  and  |i  =  cos  d,  this  can  be  written  as 

E  =  sin  a|j,  /  an  (13a) 

and  for  the  case  k'h  »  1,  that  shall  interest  us  here,  it  approaches  a  5-function  in  the  variable  p,: 

lim  sin  ap  /  ap  =  (%/a)  5(p)  (13b) 

a->Qo 

indicating  that  scattering  occurs  at,  or  close  to,  d  =  0,  or  from  Eq.  (1 1),  at  'F,=  as  stated 
previously.  From  Eq.  (10),  one  also  notes  that  the  scattering  cross  section  is  completely 
independent  of  azimuth,  as  noted  by  Ref  4.  This  indicates  that  the  scattering  takes  place  along  a 
thin  cone,  see  Fig.  4.  Its  opening  angle  is  the  same  as  the  angle  of  incidence,  and  the  scattered 
field  is  distributed  uniformly  over  the  cone.  From  Eq.  (8),  no  scattering  takes  place  if  the  incident 
wave’s  electric  field  is  perpendicular  to  the  scattering  plane;  and  likewise,  the  scattered  electric 
field  has  no  perpendicular  component.  The  physical  reason  for  these  two  facts  is  that  the  electric 
fields  need  to  have  a  component  along  the  wire  in  order  to  generate  a  current  in  it. 

For  the  case  considered  by  us,  the  width  of  the  cone,  which  is  Ap  =  7i  /  a  in  angle  for  Eq. 
(13a),  will  be  very  small.  As  a  consequence,  only  those  wires  participate  in  the  scattering  (for  the 
geometry  of  Fig.  1)  that  are  close  to  horizontal,  and  close  to  the  plane  of  incidence  and  scattering 
(the  xz  plane).  Fortunately,  the  scattering  cone  has  angle  T',  =  -T'i  so  that  the  scattered  radiation 
has  precisely  its  dominant  component  along  k  «  (Fig.  1)  where  it  is  needed. 

As  said  above,  the  scattered  field  is  purely  ^type  and  will  add  to  the  ^type  reflected  field; 
the  r-type  reflected  field  has  no  addition  from  the  scattering. 

We  can  explicitly  relate  the  amplitude  u  as  obtained  from  Huygens’  principle  to  the  near¬ 
field  scattering  amplitude  T,  using  van  de  Hulst’s  result.  Ref  7,  p.  302,  in  order  to  define  the 
phase  of  T  when  Eq.  (8)  is  used  for  this  amplitude.  In  the  latter  expression,  we  employ  'P,  = 
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as  shown  before,  and  yi  =  ys  =  0  in  order  to  obtain  the  scattering  amplitude  E|f.  Averaging  over  the 
wire  orientations  in  S  which  eliminates  the  5-function  in  Eq.  (13b)  and  furnishes  the  factor  (1/2) 
(7t/2k'h  cos  dr),  we  find 

<E'8c%  =  Eo'^t  exp(ik'r)  /{4k'r  cos  d,  [ln(0.8905  k'a)  -  i7r/2]}  (14) 

so  that  we  have  been  able  to  transform  the  cylindrically  spreading  field  of  Eq.  (8)  into  a 
spherically  spreading  field,  Eq.  (14),  as  desired.  There  still  remains  the  task  of  averaging  this 
scattered  field  over  the  azimuth  of  the  wire  orientations,  (p.  While  for  the  d-orientations,  this 
average  was  automatically  taken  care  of  by  the  5-fiinction,  it  has  to  be  done  here  explicitly  for  the 
(p  -  average. 

Fig.  5  shows  a  view  into  the  cone  opening,  indicating  its  thickness  Ad  =  7t/(2k'h  cos  d, ) 
and  since  the  cone  axis,  i.e.,  the  wire,  is  horizontal,  its  azimuthal  width  A(p.  Elementary  geometry 
furnishes  A(p  =  (2n  /  k'h)‘^ .  A  factor  2  to  be  inserted  in  the  cp-average  comes  from  the  fact  that 
wires  both  oriented  at  cp  =  0  and  9  =  7C  contribute  to  the  scattering  amplitude,  see  Fig.  1.  An 
incoherent  summation  over  all  the  N  wires  in  the  sample  provides  a  further  factor  N,  so  that 
finally  the  scattered  field  doubly  averaged  in  this  way  becomes 

<E'sc%,<p  =  E'o'exp(ik'r)  N(27c/k'h)‘'^  /{4k'r  cos  dr  •  [ln(0.8905  k'a)  -  i7t/2]}  (15) 

This  field  is  now  incident  (Fig.  1)  under  the  angle  dr  onto  the  medium  boundary,  coming  from 
inside  the  medium  and  about  to  re-emerge  into  the  air  below.  This  re-emergence  calls  for  a 
Fresnel  refraction  coefficient,  inverse  to  that  of  F^^fr  that  applied  to  penetration  from  the  air  into 
the  medium.  This  can  be  obtained  from  Eq.  (4c)  by  simply  substituting  d;  dr,  ep  ^e'p'  or 
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oranges. 

The  answer  to  this  lies  in  the  fact  that  our  sample  is  not  of  infinite  size,  but  has  a  finite 
surface  area  A.  For  a  far-distant  observer,  the  reflected  wave  fi-om  this  finite  sample  has  to  be  a 
spherically  divergent  wave  also.  This  wave  can  then  consistently  be  added  to  the  divergent  wave 
of  the  wire  scattering  process;  it  can  be  obtained  from  the  plane  reflected  wave  again  by  an 
application  of  Huygens’  principle. 

In  the  coordinate  system  (x,y,z)  of  Fig.  1,  the  reflected  wave  leaving  the  surface  equals 
Usurf=  exp  (ik"-  r)  at  z=0,  apart  from  the  amplitude  factor  E" .  The  Huygens’  formula  for  a  plane 
starting  surface  is  obtained  by  van  de  Hulst’  (p.  20)  as: 

u  =  -(ik  /  27tr)  Usurfexp(ikr)  dS,  (17) 
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and  we  assume  a  total  surface  area  of  A  =  (-dx  /  2  <  x  <  dx  /  2^  -dy  /  2  <  y  <  dy  /  2).  However,  Eq. 

A 

(17)  was  obtained  by  Ref  7  for  a  surface  normal  to  the  propagation  direction  k",  so  we  turn 
to  a  new  coordinate  system  (^,  r[,  Q  as  shown  in  Fig.  6,  and  use  the  projected  area  A  cos  Si 

as  the  starting  surface.  On  it,  the  surface  value  is  Usurf  =  exp(-ik  Q  =  1  at  Q  =  0.  Thus,  expanding 
the  radius  factor  in  the  exponential  of  Eq.  (17)  as  r  +  +  t)^)  /  2r,  we  get 


dxcos  Si/2  dy/2 

u  =  -(ik/27ir)  exp(ikr)  J  d^  j  dr|  exp  [ik(^^-hiV2r], 
-dxCOS  Si/2  -dy/2 

or 

u  =  -(2ik/7tr)  exp(ikr)  , 

with 

dxCOS  Si/2 

\  d^  exp(ik^^/2r) , 

0 

dy/2 

=  J  dri  exp(ikTi^/2r)  . 

0 

Using  the  definition  of  the  Fresnel  integrals* 

C(u)=  Jo“  cos(7tt^/2)dt, 

S(u)=  /o“  sin(7tt^/2)dt. 


(18) 

(19a) 

(19b) 

(19c) 

(20a) 

(20b) 
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we  find 


I§  =  (7tr/k)^^[C((k/7cr)  cosdi/2)  +  i  S((k/7cr)  *^dx  cosSi/2)] ,  (21a) 

I,  =  (7cr/k)''^[C((k/7tr)  '^dy/2)  +  i  S((k/7cr)  '^dy/2)] .  (21b) 

For  the  case  considered  by  us,  i.e.,  for  distant  observers,  u  is  small  and  one  has 

C(u)  =  u ,  (22a) 

S(u)  =  0(u^) ,  i.e.,  negligible, 
leading  to  the  answer 

u  =  -(ik^/27ikr)(A  cos  Si)  exp(ikr) .  (23  a) 


This  is  indeed  a  spherically  spreading  wave  as  desired,  that  can  be  added  to  Eq.  (15)  multiplied  by 
Eq.  (16);  here  Eq.  (23a)  must  be  supplemented  by  the  factor  (Fig.  1): 

Eo"'^l4iEok"'  (23b) 

TOTAL  FIELDS 

From  all  the  foregoing,  we  now  have  the  total  fields  as  seen  by  an  observer  at  large 
distance  r(carets  denoting  unit  vectors): 

(a)  perpendicular  case  (r): 

E  tot  =  E'^refl  =  Eo  E{,'  Frefl  (-ik^/27ikr)  (A  cos  di)  exp(ikr)  (24a) 

(b)  parallel  case  {C). 

Etot  =  EoEJ,"  Fjefl  (-ik^/27ckr)  (A  cos  Si)  exp(ikr) 

+  EoE,"  FJeflF''efl(27t/k'h)*^(N/4kr  cos  Sr)  [ln(0.8905  k'a)  -i7t/2]-^ 

•  exp(ikr).  (24b) 
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POLARIZATIONS 


As  stated  above,  the  correct  way  of  characterizing  polarizations  is  via  the  use  of  the  (real) 
Stokes  parameters  I,  Q,  U,  and  V.  These  are  defined  as’ 

I  =  E,E,*  +  EfEr*  =  intensity  (energy  flux/area), 

Q  =  E,E,*  -ErEr* 

U  =  E,E,*+ErE*  (25) 

V  =  i(E,E,*-E,E*) 


which  maybe  combined  into  the  Stokes  4-vector 

Q 

U 

.V. 


(26) 


If  only  linear  polarizations  are  involved,  as  is  the  case  here,  one  has  U  =  V  =  0. 
An  incident  wave  that  is  purely  linear  parallel  polarized,  has 


\.0J 


(27a) 


or  if  purely  perpendicular, 

r  n 

.0, 


(27b) 
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For  the  scattered/re-emerging  wave,  one  may  define 


1  lexpi^'^refl 

F^refr  = 

1  p'„fr  1 

exp  iT'Ufr 

(28) 

F"'^= 

1  F'U 

lexpi'F'U 

(for  a  possibly  lossy  medium),  and 

X  =  tan  [(7c/2)/ln  (0.8905  k'a)];  (29) 


further, 

I  (k"/27c)Acosdi, 

(30) 

B  «  =  I F  UF'  U I  (27c/k'h)^^  (N/4  cos  d,)[hi^(0.8905  k'a)  + 
s^„  =  X  +  ^U-Y'U  (31) 

This  yields  the  Stokes  parameters  (disregarding  the  overall  factor  1/kV): 

(a)  perpendicular  case  (r): 

Ir  =  -Qr  =  EoM  F  U  P  (k^  A  cos  (32a) 

(b)  parallel  case  (C)  \ 

=  EoM(B  refl  cos  e'refl + B  »  cos  e^sc)^  +  (B  'U  sin  e  +B  »  sin  e  K)^}(32b) 
These  equations  are  amenable  to  numerical  calculations  in  a  straightforward  manner. 
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CONCLUSIONS 


We  have  calculated  intensities  and  polarizations  of  an  electromagnetic  wave  returned  from 
a  flat,  smooth,  dielectric  sample  of  surface  area  A  that  contains  a  total  of  N  randomly  positioned 
and  randomly  oriented,  perfectly  conducting  needles.  The  returns  consist,  for  an  incident  wave 
polarized  linearly  and  perpendicularly  to  the  plane  of  incidence,  of  a  reflected  wave  only  with  the 
same  perpendicular  polarization.  For  parallel  polarization,  the  return  has  both  a  reflected  and  a 
wire-scattered  component.  The  peculiar  scattering  geometry  of  a  perfectly  conducting  wire  (in 
the  form  of  a  thin-walled  cone)  introduces  strong  reductions  in  scattering  intensity  due  to  the  fact 
that  only  few  wire  needles  satisfy  this  geometry.  However,  lesser  geometrical  restrictions  can  be 
expected  for  a  (lossy)  dielectric  wire,  with  a  corresponding  improvement  in  scattered  intensity. 

We  have  here  considered  the  case  of  wires  long  compared  to  the  wavelength,  utilizing  the 
appropriate  scattering  amplitude  of  Eq.  (8).  In  a  different  frequency  regime,  our  formalism  is 
sufficiently  flexible  to  permit  its  replacement  by  a  different  appropriate  expression  for  the 
scattering  amplitude,  including  one  corresponding  to  a  (lossy)  dielectric  wire  if  desired. 
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Figure  1.  Geometry  of  reflection,  refraction,  and  needle  scattering/re-emergence  of  waves. 


(a)  (b) 


Figure  2.  Geometry  of  incidence  on  (a),  and  scattering  from  (b)  a  thin  cylinder. 
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I 

Sc 


Figure  4.  Conical  scattering  from  a  thin,  perfectly  conducting  wire. 


AS 


Figure  5.  Vertical  and  azimuthal  widths  (AS,  Atp)  of  the  scattering  cone. 
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2 


Figure  6.  Coordinate  system  for  the  reflected  wave  spherically  divergent  fi'om  a  sample  of 
surface  area  A,  in  order  to  obtain  that  wave  using  Huygens’  principle. 
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